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ABSTRACT 

We have obtained deep Hubble Space Telescope far-UV images of 15 starburst galaxies at z ~ 1.3 
in the GOODS fields to search for escaping Lyman continuum photons. These are the deepest 
far-UV images (m^s = 28.7, 3tT, 1" diameter) over this large an area (4.83 arcmin^) and provide 
the best escape fraction constraints for any galaxy at any redshift. We do not detect any individual 
galaxies, with 3ct limits to the Lyman Continuum (~ 700 A) flux 50-149 times fainter (in f^) than 
the rest-frame UV (1500 A) continuum fluxes. Correcting for the mean IGM attenuation (factor 
~ 2), as well as an intrinsic stellar Lyman Break (factor ~ 3), these limits translate to relative escape 
fraction limits of fesc,rei < [0.03,0.21]. The stacked limit is /esc,re;(3f) < 0.02. We use a Monte 
Carlo simulation to properly account for the expected distribution of line-of-sight IGM opacities. 
When including constraints from previous surveys at z ~ 1.3 we find that, at the 95% confidence 
level, no more than 8% of star-forming galaxies at z ~ 1.3 can have relative escape fractions 
greater than 0.50. Alternatively, if the majority of galaxies have low, but non-zero, escaping Lyman 
Continuum, the escape fraction can not be more than 0.04. In light of some evidence for strong 
Lyman Continuum emission from UV-faint regions of LBGs at z ~ 3, we also stack sub-regions of 
our galaxies with different surface brightnesses, and detect no signficant Lyman continuum flux at 
the fesc,rei < 0.03 Icvcl. Both the stacked limits, and the limits from the Monte Carlo simulation 
suggest that the average ionizing emissivity (relative to non-ionizing UV emissivity) at z ^ 1.3 is 
significantly lower than has been observed in Lyman Break Galaxies (LBGs) at z 3. If the ionizing 
emissivity of star-forming galaxies is in fact increasing with redshift, it would help to explain the high 
photoionization rates seen in the IGM at z > 4 and reionization of the intergalactic medium at z > 6. 
Finally, we see no evidence of shadowing of the far-UV background by Hi in our target galaxies 
(on scales smaller than ^ 40 kpc), suggesting that either the majority of the far-UV background 
originates at lower redshift (or within our galaxy) or that the high column-density Hi in our target 
galaxies is distributed over larger than 40 kpc scales. 

Subject headings: galaxies: high- redshift, galaxies: starburst, galaxies: intergalactic medium, ultravi- 
olet: galaxies 



1. INTRODUCTION 

The Hi-ionizing background and its evolution are crit- 
ically important to various aspects of galaxy assem- 
bly: reionization of the intergalactic medium (IGM) by 
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z ~ 6 (|Fan et al.ll200l : supression of star formation in 
low mass dark rnatter halos (< 10^ Mm') (iRees 119861 : 
lEfstathioul [l99l IBuUock et all \20M. ISomervilld 12001) : 
and affecting the Hi distribution in outskirts of galaxies 
povc & Shuiri994). Unfortunately, the sources of the 
ionizing background at high redshift are not yet known. 
Searches for high redshift (z > 3) QSOs show that their 
space densities fall precipitously toward higher redshift 
and are too rare at the highest redshifts t o contribute sig- 
nificantly t o the Hl- ionizing background (Richards et aLl 
2006; Ho pkins et ah 2007; Siana ct al. 2008; Jiang et al] 
20081 i2009t ). Therefore, young star- forming galaxies have 
become the most likely candidates for providing the ion- 
izing emission. However, it isn't clear what fraction of 
the ionizing photons (hereafter called the escape fraction, 
fesc) produced by massive stars can escape the high col- 
umn density Hi gas surrounding these star-forming re- 
gions. Many studies have been conducted searching for 
escaping ionizing radiation from galaxies at various red- 
shifts, with most providing only uppe r limits to the ion- 
izing emissivity. iSteidel et al.l (|2001[ ) reported the first 
detection of Lyman continuum (LyC) photons in a stack 
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of 29 z '-^ 3.4 LBG spectra. The flux ratio between the 
rest- frame 1500 A and the Lyman continuum (corrected 
for average IGM absorption), /1500//900, imphes that the 
majority of Lyman Continuum photons produced by the 
massive stars are escaping into the IGM. Deeper spectra 
of 14 LBGs show that most have low relative escape frac- 
tions {fesc,rei < 0.25), with & Sma ll fraction (2/14) ex - 
hibiting esca pe fractions near un ity (jShaplev et al.|[2006l ). 
In addition, llwata et"an (|2009f ) have conducted narrow- 
band imaging of the Lyman continuum of 125 Lya emit- 
ters (LAEs) and 73 Lyman break galaxies (LBGs) at 
z = 3.09 and detect ~ 10% with large escape fractions. 

While many galaxies at 2 ~ 3 have now been 
detected in the Lyman continuum, all searches at 
low red shift (z ~ a nd z ^ 1) have g iven null 
results (iLeitherer et al.l [19951 iDeharveng e t al. 2001j 



Giallongo et al. 2002; Malkan et al.l '2003; 



2007tlGrimes et al.l2007:.Cowie et al.t2009„: Grimes et al.l 



2009( 1. The most significant results have come from 
sgace-based far-UV imaging with HST and GALEX. 
Siana et al. (2007) compiled (along with data from 
Malkan et al.l I2003il a sample of 32 galaxies at 1.2 < 



z < 1.5 with deep HST far-UV images with no sig- 
nificant detection of escaping Lyman continuum. This 
implies that the galaxies must, on average, have a UV- 
to-LyC continuum ratio o f /i5oo//Li/c(3cr) > 50. Re- 
cently, (jCowie et al.ll2009 t) stacked GALEX far-UV im- 
ages of 626 star-forming galaxies at 0.9 < z < 1.4 
and did not get a signficant detection, with a limit of 
/i50o//Lyc(3o-) > 83. These studies suggest that ei- 
ther the average escape fraction has been decreasing with 
time, or that analogous galaxies at low redshift have not 
yet been targeted. However, given the small fraction of 
galaxies at z ~ 3 with observed large escape fractions 
(~1/10), we still would only expect a few detections thus 
far at z ~ 1 if their escape fractions are similar to higher 
redshift galaxies. Thus, larger samples are needed at low 
redshift. 

Most published escape fraction limits at z 3 are not 
very deep, with typical /esc, re; (^c) 0.25. Therefore, it 
is possible that many of the galaxies with non-detections 
do emit significant ionizing continuum below current de- 
tection limits. Given that there are 10 times more ob- 
jects with non-detections, even small amounts of ioniz- 
ing emission from these galaxies would contribute signif- 
icantly to (or even dominate) the ionizing background. 
Therefore, deep surveys probing low escape fractions are 
also important in determing the total contribution of star 
formation to the ionizing background. 

In addition to simply measuring the escape fraction 
and its evolution, it is important to understand how 
these photons are able to escape into the IGM, given 
the presumably large reservoirs of gas and dust in these 
actively star-forming galaxies. Several mechanisms have 
been invoked to explain high escape fractions, each of 
which predicts unique observational signatures. Both 
semi-analytic models and numerical simulations suggest 
that SNe winds may be sufficiently powerful to expel gas 
from the disk, producing low Hi column density lines 
of sight (LOSs), though it is highly dep endent on the 
spati al distribution of the star clusters ([Clarke fc OevI 
I2OO2I ) an d the star formati on efficiency (relative to the 
gas mass. [Fuiita et al.[l2003f ). In such a scenario, escaping 



Lyman continuum may be expected in only the regions 
of the galaxy with high star formation surface density . 
The simulations of iRazoumov fc Sommer-LarsenI (|2006l 
[S)07, 2009) suggest that fe edback is far n i ore ef ficient in 
lower mass halos. However. I Gnedin et al.l (2008) suggest 
that the smaller gas-to-stellar scale height ratios in high 
mass disks results in a higher probability o f young stars 
being exposed to low column LOSs (see also [Wise &: CenI 
[20091 for high redshift dwarf galaxy results). In this sce- 
nario, we would expect to see escaping LyC photons from 
above (in azimuth) an ed ge on disk and only in the most 
massive galaxies. Finallv. TRicottil ([2002[ ) posits that it is 
the formation of globular clusters (GCs) that may have 
reionized the universe. In which case, we would expect 
to see no evidence for large escape fractions at z 1, as 
the vast majority of GCs had been formed long before 
then and should no longer contain massive stars. 

Determining from which parts of the galaxies the pho- 
tons are escaping can help explain the mechanisms by 
which such large fractions of Lyman continuum photons 
can escape these large Hi reservoirs. Does the ionizing 
emission simply trace the UV (1500 A) emission, or is it 
only escaping from certain portions of the galaxy? What 
are the spectral energy distributions (SEDs) of the LyC 
emitting regions? And if the LyC is only escaping from a 
portion of the galaxy, how can such large ionizing fluxes 
observed in some LBGs be emitted from only small frac- 
tions of the total starburst? 

With these questions in mind we conducted a deep, 
high spatial resolution, far-UV imaging program of the 
15 most UV- luminous star-forming galaxies at 1.2 < 
z < 1.5 in th e Great Obser vatories Origins Deep Sur- 
vey fGOODS. [Giavalisco et al. 2004). We use the Solar 
Blind Channel (SBC) of the Advanced Camera for Sur- 
veys (ACS) on the Hubble Space Telescope (HST) with 
the F150LP filter, which samples the rest-frame Lyman 
Continuum ('^ 700 A) at these redshifts. These data are 
typically a factor o f two deeper than our previous survey 
([Siana et al.l |2007[ ) and target galaxies about a magni- 
tude brighter (at Xrest ~ 1500 A), giving a factor of 5-10 
higher sensitivity to escaping LyC. Also, the high spa- 
tial resolution allows us to investigate variations in the 
escape fraction via measurement of the relative distribu- 
tion of the Lyman continuum and rest-frame UV fluxes 
to help determine the mechanism by which the Lyman 
continuum photons are allowed to escape. 

Throughout the text we use a flat ACDM cosmology 
with Ho = 70 km s'^ Mpc-\ Qm = 0.3, and = 0.7. 
All flux densities are given in f^, and all magnitudes are 
in the AB system. 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. Target Selection 

We targeted the 15 most UV- luminous galaxies in 
the GOODS fields with spectroscopic redshifts between 
1.2 < z < 1.5. The sample has the brightest GOODS 
B-band (F435W, Xrest ^ 1900 A) fluxes as weU as the 
bluest B — I colors, in order to avoid dusty and/or 
older galaxies. The specified redshift range keeps the 
Lyman limit redward of any significant system through- 
put (z > 1.2) and ensures high sensitivity to the rest- 
frame UV continuum (z < 1.5). These galaxies have 
luminosities typical of the L^y LBGs at z ~ 3 that ex- 
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hibit large escape fractions. The targets were chosen in 
the GOODS fields for the wealth of ancillary data (HST, 
Spitzer, Chandra, optical spectroscopy, etc.) necessary 
for analyzing various galaxy properties (star formation 
histories, stellar masses, metallicities) . Finally, all X-ray 
sources were excluded to avoid contamination by unob- 
scured AGN. 

2.2. Observations 

We imaged each galaxy with the Solar Blind Channel 
(SBC) of the Advanced Camera for Surveys (ACS) and 
the F150LP filter, a blue cutoff filter blocking light with 
A < 1450 A. This setup effectively blocks the geocoronal 
Lya and [Oi] emission that would significantly increase 
the background in our images. The red cutoff of the 
throughput is dictated by the decreasing sensitivity of 
the Multi- Anode MicroChannel Array (MAMA) towards 
redder wavelengths. The configuration gives very low 
throughput at A > 2000 A (< 0.0008 of peak sensitivity 
at 1500 A), with further decline in sensitivity at red- 
der wavelengths. We note here that the rate of decline 
in the QE of the MAMA at A > 2000 A has recently 
been found to be shallower than previously thought, re- 
sulting in a higher than expected QE at optical wave- 
lengths. However, we have performed tests with the new 
throughput curve for various SED scenarios and find that 
leakage of light redward of A > 2000 A is not a con- 
cern in this study as the expected leaking flux is below 
the noise limits achieved in this study. The transmis- 
sion weighted effective wavelength is Ae// = 1610 A, or 
Arest 700Aat z ^ 1.3. We note that because the pho- 
toionization cross-section falls as approximately v^^ , the 
HI opacity at Xrest ^ 700 A is nearly half of that at 900 

A. 

Each galaxy was imaged for five orbits, split into two 
visits of two and three orbits to ensure that SBC was not 
in use for too long, as th e dark current incre ases as the in- 
strument warms up fsee lTeplitz et al.|[2006[) . In addition, 
we also ensured that the visits did not follow any previ- 
ous SBC use to ensure a low MAMA temperature. Each 
orbit consisted of four dithered exposures of 630 seconds, 
giving a total exposure time of 12600 seconds per galaxy. 
Finally, because the dark current glow is concentrated 
near the center of the detector, we placed the objects in 
the lower-left quadrant of the detector where the dark 
current glow is at a minimum. 

2.3. Data Reduction 

First we masked all known bad pixels and then sub- 
tracted the most recent dark image provided by Space 
Telescope Science Institute (STScI). This dark image 
does not reflect the additional (and often much stronger) 
dark current "glow" that arises when the MAMA warms 
up. We use our data to produce a profile for this dark 
"glow" to be subtracted from each image. To do this 
we mask all known objects in the GOODS B-band im- 
ages usmg_tlie_segmcntation maps output from SExtrac- 
tor (|Bertin fc Arnou ts 1996) and add all of the exposures 
together. A two dimensional 4th order spline is fit to the 
stacked dark, giving a smooth profile for the dark "glow" . 
The amplitude of the "glow" changes with each exposure, 
so we scale the dark glow fit to the background of each 



image (with known objects masked) and subtract it. The 
final, dark subtracted frames are then fiattened using 
the fiat field provided by STScI. Th ese reduced frames 
are drizzled (|Fruchter fc Hookl [200l ) onto the GOODS 
tiles to facilitate direct comparison with the other data 
{pixfrac — 0.6, output pixel scale= .03"). The relative 
astrometry between objects detected in both the SBC 
and GOODS B-band (F435W) images is less than one 
pixel (0.03"). 

In the region of the images where our primary targets 
are located, there is very little of the additional dark 
current glow, and the resulting drizzled images reach a 
3cr sensitivity of uiab = 28.7 in a 1" diameter aperture. 
The regions near the center of the detector can be up to 
one magnitude less sensitive due to the additional dark 
current (which dominates the noise). These data are the 
deepest images ever obtained at these wavelengths and 
over this large an area (15 x 0.32 — 4.83 sq. arcmin). 

To define apertures for far- UV flux determinations , 
we started with the SExtractor (jBertin fc Arnoutill996[) 
segmentation maps of the GOODS z' (F850LP) images. 
These apertures are much larger than the areas with 
significant emission at rest-frame UV wavelengths. We 
therefore chose to shrink the apertures to isophotes that 
contain 90% of the GOODS B-band (F435W). This re- 
duces the aperture sizes by about 50%, resulting in sen- 
sitivities about 0.4 mags deeper. Using the B-band data 
to determine the aperture is reasonable, as it is unlikely 
that a source emitting strongly at Xrest 700 A is not 
also emitting strongly at Xrest ~ 2000 A. 

Finally, we have corrected the far-UV measurements 
for Galactic extinction using the estimates from the 
100 fim COBE/DIRBE maps {Ay =0.0 4 and 0.026 in 
GOODS-N and GOODS-S resoe ctivelv. ISchleeel et al.l 
1 1998. ). The extinction curve of ICardelh et al.l ()1989f ) 
with an Ry = 3.1 gives an AiQio/Ay = 2.55, result- 
ing in ^1610 =0.10 and 0.07 mags, respectively. As these 
fields are some of the lowest extinction site lines out of 
the Galaxy, the extinction corrections are less than 10%, 
even in the far-UV, and any uncertainties in the exact 
shape and normalization of the extinction curve do not 
significantly affect the photometry. 

3. RESULTS 

3.1. Individual Limits 

The final far-UV images of the galaxies are shown in 
Figure m alongside the GOODS B-band (F435W) images. 
None of the ta rget galaxies were detected in the SBC 
far-UV images. Ilwata et"aLl ()2009f ) see significant offsets 
between the detected LyC and the rest-frame UV con- 
tinuum (typical offset ~ 1", or ^ 7.6 kpc at z — 3.1). 
Therefore, we have also searched for escaping Lyman 
continuum within a 1.5" radius of our targets. No de- 
tections are seen by eye. In Table [T] we list the 3cr limits 
to the Lyman continuum flux, fhyC^ and the UV-to- 
LyC continuum flux ratio limits, f 1500/ fLyC- The 3a 
limits to the UV-to-LyC flux ratios (in /^) are shown 
in Figure [5] and range from 50-174. The depths of the 
images are nearly uniform, so the large range in UV- 
to-LyC limits is caused by the factor of 3-4 variations 
in both aperture area and 1500 A flux. Before compar- 
ing these /i5oo//lj/C limits to other surveys, it is impor- 
tant to correct for the average IGM opacity (factor of 
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Fig. 1.— GOODS B-band (F435W, left) and far-UV (F150LP, right) images of all 15 targets. The stamps are 4" on a side (1" = 8.4 kpc 
at 2 = 1.3). The target galaxies display a variety of morphologies: including spirals, compact galaxies, mergers, and/or clumpy disks. 



- 2, iSiana et alJ 120071 : llnoue fc Iwatall2008f ) as this is a 
strong function of both redshift and wavelength. These 
IGM-corrected hmits, {fLyc/fi50o)cor, can now be di- 
rectly compared to other IGM-corrected detections and 
limits of other surveys. We choose to use this ratio for 
comparison, rather than an estimate of the escape frac- 
tion, as there are additional assumptions made when de- 
termining the escape fraction. However, for comparison 
with other surveys, we will also convert this to a "rela- 
tive" escape fraction, which is defined as the fraction of 
LyC photons escaping the galaxy divided by the fraction 
of escaping 1500 A photons. The relative escape fraction 
is computed as 




fe 



(/l50o/ Jlc) 



LCjstel 



if 150o/ fLc)obs 



exp[TiGM) 



(1) 



where tjgm ^ 0.6 is the LyC optical depth of the IGM 
at z = 1.3 (as measured through the F150LP filter, 
ISiana et all 120071) and (/i5oo//lc) stei = 3 is the intrin- 
sic Lym an Continuum brea k of the stellar population (in 
/^). In ISiana et all ((20071) we argue that given typical 
metallicities and star formation histories for these galax- 
ies, the intrinsic break may be much larger (6 — 10) but 
use this value to be consistent with the values assumed 



50 100 150 0.0 0.10 0.20 0.30 
fi500 / fiyc Limits (3ct) f„,_,,| Limits (3<t) 

Fig. 2. — Histogram of the 3cr limits to the UV-to-LyC ratio, 
f 1500/ f LyC > (left) and the resulting upper limits to the ionizing 
escape fraction (right). To determine the relative escape fraction 
limits an intrinsic Lyman break of 3 (in f^) was assumed and a 
redshift dependent IGM correction (factor of ~ 2) was applied. 
The arrows indicate the direction of the limits, and emphasize that 
these are not detections. 

in other works. 



3.2. Stacked Image 
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TABLE 1 

Target Information 



Name 


RA (J2000) 


Dec (J2000) 


^spec 


U 


Area'' 






fl500/ fLyC'^ 


f 

J esc,rel 


s4 


03:32:11.57 


-27:46:22.9 


1.221 


23.48 


3.64 


1.66 


<0.0182 


>90 


<0.06 


s7 


03:32:13.66 


-27:43:13.1 


1.856 


23.46 


1.50 


1.64 


<0.0138 


>118 


<0.21 


sll 


03:32:23.96 


-27:43:49.1 


1.311 


23.32 


1.69 


1.59 


<0.0091 


>174 


<0.03 


si 7 


03:32:37.37 


-27:50:13.7 


1.389 


23.26 


2.06 


1.85 


<0.0144 


>128 


<0.05 


s23 


03:32:46.79 


-27:52:29.8 


1.227 


23.56 


2.82 


1.26 


<0.0180 


>69 


<0.08 


slO 


03:32:22.92 


-27:48:57.4 


1.264 


23.47 


4.02 


1.41 


<0.0208 


>67 


<0.08 


s6 


03:32:13.21 


-27:41:58.0 


1.297 


23.73 


2.67 


1.12 


<0.0182 


>61 


<0.09 


sl5 


03:32:35.25 


-27:54:32.6 


1.414 


23.66 


1.22 


1.22 


<0.0082 


>149 


<0.05 


s8 


03:32:18.54 


-27:48:34.1 


1.414 


23.75 


3.81 


1.00 


<0.0198 


>50 


<0.14 


nl 


12:36:12.42 


62:14:38.4 


1.432 


22.90 


5.17 


2.96 


<0.0200 


>148 


<0.05 


ii3 


12:36:28.94 


62:06:16.0 


1.265 


23.25 


7.13 


1.34 


<0.0216 


>61 


<0.09 


n5 


12:35:55.87 


62:13:32.7 


1.296 


23.25 


3.12 


1.52 


<0.0164 


>92 


<0.06 


nil 


12:37:05.59 


62:17:17.8 


1.250 


23.60 


1.95 


1.24 


<0.0125 


>98 


<0.06 


n6 


12:36:45.71 


62:07:54.3 


1.433 


23.92 


4.09 


1.04 


<0.0207 


>50 


<0.14 


nl4 


12:37:09.12 


62:11:28.6 


1.341 


23.63 


1.92 


1.19 


<0.0148 


>80 


<0.08 


Stacks 
High SB <= 
Low SB ° 












9.31 
11.5 


0.0199 
0.0575 


> 469 

> 200 


<0.013 
<0.030 



arcmin 

0.9 X the total value (to properly match far-UV aperture) 
Limits are 3cr 

SB — surfaee brightness, high (low) SB stacked above (below) Sb{AB) — 23.06 mag arcsec" 



We have stacked both the F435W and FUV images of 
the 15 targets with a simple addition of all of the images 
(See Figure 131) ■ No obvious detection is seen in the stack. 
Of course, due to the varying sizes and morphologies of 
these galaxies, some galaxies are adding noise to areas 
where other galaxies have strong B-band flux. Therefore, 
we do an optimized stacking which consists of adding up 
only the pixels that were used in the individual aperture 
extractions (based on the GOODS F435W isophotes). 
We sum all of the pixels in the individual galaxy aper- 
tures (and their associated errors in quadrature) and get 
a total LyC flux of fLyC = 2.2 x IQ-^i ± 2.0 x lO^^i 

cm~^ Hz~^, not a significant detection. The limit 
to the UV-to-LyC ratio is /isoo/ZiyC > 302 (3ct). After 
correcting for the average IGM transmission, this gives a 
3(T limit to the LyC flux 151 times fainter than the UV 

flux^ 

iCowie et al.l ([2009I ). in their GALEX far-UV stack of 
626 galaxies, got a ~ 2cr negative flux in the location 
of the stack relative to the background. Though it was 
unclear whether this may have been a systematic effect, 
it was suggested that a negative flux could be caused by 
a shadowing of the far-UV background by the Hi in the 
target galaxies. In our stacked image we see no evidence 
for a flux below the background level. In Figure 2] we 
show a zoomed out version of our far-UV stack smoothed 
with a 0.5" Gaussian. There is no evidence for a shadow 
in this image ( with scales less than 5", or 40 kpc). The 
shadow seen bv ICowie et al.l (|2009| ) appears to be large 
(5-10"), but the signal is very low, and this is of order 
the size of the GALEX far-UV beam so it is unclear 
whether the absorption is real or over what angular scale 
the absorption is taking place. The lack of shadow seen 
in our data i s perhaps not surpr ising, as most of the UV 
background ("Brow n et al.l l2000l ) likely originates in the 
foreground of these g alaxies, either from Galactic sources 



(jHurwitz et al.l 



galaxies ()Henrvlll99lD . 



19911) or Lya emission from foreground 



F435W Stack 



Far-UV Stack 




It is still not known which parts of galaxies are likely to 



Fig. 3.— The B-band (F435W, left) and far-UV (F150LP, right) 
stacked images. The images are a simple sum of the individual 
objects. No detection is seen in the far-UV image. 



have high LyC escape fractions, or what physical proper- 
ties dictate the escape fraction. One property that may 
affect the escape fraction is the star-formation surface 
density. On one hand, regions of high star-formation sur- 
face density may provide sufficient winds to expel gas and 
produce low Hi column density lines-of-sight. Alterna- 
tively, regions of low star formation surface density (like 
the far outskirts of disks in the local universe) have lower 
Hi columns so feedback may therefore be m ore effective 
at dis rupting the gas distribution. Indeed, llwata et"al] 
()2009() suggest that many of their LyC detections are 
offset by several kiloparsecs from the brightest regions of 
the galaxies. It is unclear whether this emission is coming 
directly from star-formation or from nebular free-bound 
continuum emission fro m ionized gas on the outskirts of 
the galaxy ()Inouell20"09l ) . To determine if such extended 
LyC emission exists in our galaxies, we first checked all of 
the individual images by eye to make sure that our aper- 
tures (defined by the rest-frame 1900 A continuum) were 
not missing any LyC flux with extreme offsets (> 1"). 
Second, within our previously defined B-band apertures, 
we chose to separate regions with low and high UV sur- 
face brightness (dividing the sample at Sb{AB) = 23.06 
mags arcsec"^, the typical lOc sensitivity of a single 
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Fig. 4. — Same far-UV (F150LP) stacked image as in Figure 
[S] but zoomed out by a factor of two and smoothed with a 0.5" 
Gaussian. No clear detection exists at the center, nor is there 
a "shadow" of the background by our target galaxies. Detected 
galaxies have not been subtracted in order to give an idea of the 
low level of confusion. 

pixel in the GOODS B-band images (corresponding to a 
star formation surface density of 0.187 Mq yr~^ kpc~^). 
Both the high and low surface brightness regions give 
non-detections in the LyC. {fisoo/ fLyc)high-sb > 469 
(3cr) and {f 1500/ fLyc)iow-sb > 200 (3cr). Note that we 
do not have resolved U-band images of these galaxies so, 
for these surface brightness stacks, we used the B-band 
fluxes and multiplied by 0.85 (the average flux ratio be- 
tween the U and B images) to estimate /isoo- The details 
of these stacked limits are also provided in Table [TJ 

4. COMPARISON WITH STUDIES AT Z ~ 3 

In Figure [5] we show the f LyC / f 1500 measurements 
(versus UV luminosity) from several large surveys at 
z ~ 1 and 2: ~ 3 corrected for the average IGM at- 
tenuation at the relevant wavelengths (Ai-est ~ 700 A at 
z ^ 1.3 and Xrest 900 A at 2; 3). 

At 2: ~ 3, recent studies have shown large LyC- 
to-UV ratios ( f^,vc/ f^F,nn > 0.2) in - ^10% of galaxies 
(jShaplev et al.l 120061 [l^ata et all 120091) . However, there 
are no LyC detections at z 1 despite the large number 
of galaxies (47) that have 3a upper limits to /lj/c//i500 
that are lower than any ratios seen in z > 3 galaxies with 
LyC detections. 

It appears evident that a higher fraction of galaxies 
have large escape fractions at z ~ 3 than at z ^ 1.3. 
However, it is first important to determine if these galax- 
ies have similar stellar populations, as differences in the 
stellar populations will affect the intrinsic LyC fluxes 
(and the implied escape fractions). Also, the charac- 
teristics of the foreground IGM are different in the two 
samples for two reasons: the number of absorbers (per Hi 
column density) drops dramatically between z ^ 3 and 
z 1 and different rest-frame wavelengths are probed at 
the two redshifts (700 A at z ~ 1.3 and 900 A at z ~ 3). 
Therefore, in order to properly compare the samples at 
z ~ 3 and z ~ 1, we first need to determine that the 
two galaxy populations have similar stellar populations 



(and thus intrinsic /lj/c//i500 ratios), and then deter- 
mine the effects of the IGM through Monte Carlo simu- 
lations based on the known redshift and column density 
distributions of Hi absorbers in the IGM. 

4.1. Galaxy Characterization 

To derive an escape fraction from the /iyc/Zisoo 
limits, we must first have an estimate of the intrin- 
sic f LyC / f 1500 ratios. Specifically, we must be certain 
that the stellar populations of these galaxies have strong 
intrinsic Lyman Continuum emissionand have similar 
properties to the z ~ 3 LBGs that have been shown to 
exhibit high e scape fractions. To these ends, we fit the 
iBruzuall (|2007l ) stellar population models to the available 
optical (ground-based U, GOODS HST B,V,I,z), near- 
IR, and Spitzer 3.6 and 4.5 /im photometry to estimate 
the intrinsic LyC flux, stellar mass, starburst age, and 
dust reddening for each galaxy. 

For each galaxy we assume a Sa lpter IMF (iSalpeteii 
Il955| l. a Calzetti reddening curve (ICalzettil[l997r and 
solar metallicity. We fit to three different exponentially 
declining star formation histories (SFR oc e"*/"^) with e- 
folding times of r =0.1 and 0.3, and oo (constant SFR) 
Gyr. The ages, star formation rates, and dust reddening 
are allowed to vary. 

Though the S/N of the photometry is high and spans 
the rest-frame UV, optical and near-IR, the precise star- 
formation history is still ambiguous. That is, the val- 
ues are similar for the best-fit SEDs to the three differ- 
ent star-formation histories. This is not too surprising, 
as the effect of the subtle differences in star-formation 
history can be masked by varying reddening and age. 
Unfortunately, the best-fit SEDs using the three differ- 
ent SFHs give very different predictions of the intrin- 
sic Lyman Continuum fiux. This is because the portion 
of the spectrum to which we are fitting our photome- 
try is dominated by stars with ages > 10 Myr (typically 
greater than 100 Myr), whereas the Lyman Continuum 
flux is dominated by massive stars that have lifetimes 
shorter than 10 Myrs. This is demonstrated in Figure 
[6l where the flts with three different assumed SFHs are 
plotted. At wavelengths where we have photometry, the 
difference between the fits is negligible. However the in- 
trinsic SEDs (ie. not attenuated by dust or Hi in the 
IGM, shown as dotted lines) exhibit very different Ly- 
man Continuum fluxes (with typical variatons of a fac- 
tor of ^ 5 for our galaxies). However, we think it is 
important to note that a constant star formation his- 
tory gives a maximum Lyman Continuum flux of any of 
the typically considered histories (single burst, constant, 
or exponentially decline) and using this SFH still yields 
Lyman-break amplitudes of /isoo/.f^f 8 , similar to 
what was found in iSiana et al l (120071) . I n Sian a et al.l 
(j2007D we fit to both thelBruzual fc Ch"arlot (200§} and 
Starburst99 models (jLeitherer et al. 1999) and got very 
similar Lyman-break amplitudes and stellar population 
parameters. 

Regardless of how well we can predict the intrinsic LyC 
luminosities of these galaxies, the SED fitting is useful for 
comparing our sample to higher redshift galaxies selected 
via the Lyman Break technique. In Figure [71 we compare 
ou r best-fit parameters to those of a sample of 72 LBGs 
of iShaplev et al.l (|2005f) . In both cases, a constant star 
formation history, Calzetti reddening law, Salpter IMF, 
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Fig. 5. — Our new fLyc/fl500 measurements with other individual limits at 2: ~ 1 l|Malkan et al.|[200 3': 'Siana ct al. '2007") and z ~ 3 
llShapIev et al. 2006; Iwata et ani2009l) . The vertical dashed and dotted lines represent M*^qq at z ~ 1 (-20.2, Arnouts et al. 2005), and 
2 ~ 3 (-20.97. IReddv et al.l l200S). respectively. The studies at both redshifts sample similar absolute magnitudes. All reported limits have 
been adjusted to 3cr and each limit has been corrected to account for the mean IGM transmission at each redshift. Nearly all of the z ~ 1 
measurements have f LyC / f 1500 limits below the observed ratios of z ~ 3 galaxies with LyC detections. 



and solar metallicity were assumed. The distributions of 
dust reddening, SFR, and stellar mass are very similar to 
the distributions for the z ~ 2 LBGs, though the best- 
fit stellar ages span a narrower range toward the low 
end of the LBG distribution (at ~ 300 Myr). Given the 
similarities in the derived parameters, we are confident 
that any differences in the LyC-to-UV flux ratios between 
these two samples are not due to fundamental differences 
in the stellar populations. 

4.2. Monte Carlo IGM Modeling 

In Figure [HI we corrected for the average IGM attenua- 
tion of the LyC at the relevant wavelengths. At z ~ 1.3, 
the mean IGM transmission is about 0.5 bu t the distri- 
bution along different line s-of-sight is large (jSiana et al.l 
[2007LlInoue &: Iwatall2008D . In 10-15% of LOSs, the IGM 
is completely opaque and we would not expect to see any 
flux, regardless of the escape fraction. Alternatively, 25- 
30% of LOSs have transmission greater than 0.8 and our 
limits are far more sensitive than assumed in Figure [S] In 
order to properly consider the large range in IGM optical 
depths towards galaxies at z 1.3, we have performed a 
Monte Carlo simulation to determine the escape fraction 
parameter space in which our galaxies must lie in order 
to get non-detections for every object. Th e method is 
exactly the same as the method outlined in iSiana et al.l 
(U007), but is briefly summarized here. 

In our simulations, we assume that all of these blue, 
UV-luminous galaxies have an intrinsic flux ratio (in /i,) 
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Fig. 6. — Stellar population model fits to the optical, near-IR, 
and mid-IR photometry. Three star formation histories are used 
(constant star formation and exponentially declining star forma- 
tion with e-folding times of 100 & 300 Myr. The amplitude of the 
Balmcr break is very small, indicating young < 300 Myr old star- 
bursts. Note that, although the fits with different SEDs give very 
similar values, the predicted intrinsic Lyman Continuum varies 
significanty (typically by a factor of ~ 5 in our fits). 

of /i5oo//lj/C = 3. We assume a simplified parameter 
space for the relative escape fraction such that X frac- 
tion of galaxies have a relative escape fraction, Y, and 
the rest have fesc.rei = 0. We step in X,Y parame- 
ter space in 50 evenly spaced intervals (between and 
1). For each location in this parameter space, we ran- 
domly select 10,000 galaxies from this assumed distri- 
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Fig. 7. — Comparison of galaxy param eters for our z ^ 1.3 s ample (thick line) compared with the sample of l lShaplev et al.ir2005l . shaded 
histogram). For ease of comparison, the 'Shapl ev et all 120051) histograms have been scaled so that they have the same normalization as 
our sample (actual sample contaions 72 galaxies). Though our sample is small, it can bee seen that our galaxy sample has very similar 
properties to the higher redshift (z ~ 2) LBG samples. 

bution and place them along lines-of-sight with opaci- 
ties rand omly selected from the IGM distributions deter- 
mined in iSiana et al.l ()2007[ ) . We then determine what 
fraction of the time our experiment would produce only 
non-detections (as we observe). The parameter space 
where we would obtain non-detections 84.1, 97.7, and 
99.9% of the time are considered our 1, 2, and 3a lim- 
its on the escape fraction parameter space. We use the 
list of all the observed /isoq/Z lmC 3(7 limits from ou r ob- 
servations, as_well as those of iMalkan et al.l ()2003D and 
ISiana et al.l (|2007D . 

The results of this Monte Carlo simulation are summa- 
rized in Figure jS) As can be seen, the vast majority of 
the parameter space is ruled out with our observations. 
At the 95% confidence level there can be no more than 
8% of galaxies with fesc,rei > 0.5 or no more than 50% of 
galaxies with fesc,rei > 0.04. Both the increased depth 
of the survey and the larger sample help sig nificantly in- 
crease the constraints over those achieved in ISiana et"al] 
(|2m7[) . 

Unfortunately, a similar analysis is difficult for LBGs 
and Lya-emitters at z ~ 3 . First, ther e is a c oncern that 
some of the detections by llwata et al.l (|2009[ ) may actu- 
ally be due to contamination from foreground sources. 
Some of these LyC detections have large offsets (some 
greater than 10 kpc) from the UV-luminous LBGs with 
spectroscopic redshifts, raising concerns that the pre- 
sumed LyC emission is actually non-ionizing emission 
from faint foreground galaxies. 



In addit ion, many o f the L yC detected Lya emitters 
(LAEs) in llwata et all (|2009f ) display LyC fluxes (in f^) 
equal to or greater than the rest-frame UV (1500 A) con- 
tinuum, whereas stellar population models display a large 
decrement shortward of the Lyman limit. Ilnoud ()2009f ) 
argue that this LyC "bump" may be due to free-bound 
continuum emission from the surface of dense Hi regions 
that are illuminated by star-forming regions with large 
escape fractions. This mechanism would result in more 
energetic LyC photons being reprocessed by the Hi re- 
gions and re-emitted (via free-bound emission) at wave- 
lengths just short of the Lyman Break. If the electron 
temperature of the gas is fairly low 10'* K), the ther- 
mal energy of the electron is insufficient to add signifi- 
cant energy to the resulting bound-free photons. Thus, 
this LyC "bump" would not extend to wavelengths much 
below the Lyman limit 750A). This would naturally 
explain why we do not see such strong LyC fiux with our 
observations, because we are most sensitive to LyC pho- 
tons at lower rest-frame wavelengths (650-750 A). How- 
ever, we do not believe that free-bound nebular emis- 
sion is significant in our sample at z 1.3, because 
rather unique conditions are required (extremely young 
ages (< 1 Myr), low metallicities (<l/50 Zq), and/or 
top-heavy IMFs), none of which should be prevalent in 
massive starbursts at z ~ 1.3. If some of the z ~ 3 LBG 
LyC detections are actually contaminants or if the 900 
A LyC flux in LAEs is boosted via free-bound emission, 
this makes direct comparison with our results difficult. 
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Fig. 8. — Escape Fraction parameter space at 2 ~ 1.3 allowed by 
our o b servations (i n com bination with limits from IMalkan et al.l 
II2003I ): [Siana et al.l 120071) ). In our simplified parameter space, a 
fraction of galaxies (y-axis) have a non-zero relative escape frac- 
tion (x-axis). The darker shaded areas show the parameter space 
excluded at the 1, 2, and Scr levels from our Monte Carlo simu- 
lations. At 95% confidence level, we can rule out 8% of galaxies 
having /^jj. ^.^i > 0.5 or, alternatively, rule out 50% of galaxies 
having fesc ',rel < 0.04. 

as it is unclear what the intrinsic LyC flux is in either 
case. 

Generally, it appears that z ^ 3 studies are finding 8- 
10% of g alaxie s wit h large escap e fractions iShaplev et alj 
(1/14 in [20061) and llwata et"al] (7/73 LBGs and 10/125 
LAEs in 2009). This is also the case in a much larger 
spectroscopic survey of z 3 LBGs, where ^^10% are 
found with large escape fractions (M. Bogosavljevic 2009, 
private communication) . The results plotted in Figure [S] 
show that, if 8% oi z ^ 1.3 galaxies have large relative 
escape fractions (> 0.5), then we should have seen at 
least one detection in the z ~ 1.3 studies in 95% of our 
IGM realizations. 

Finally, we note that in iSiana et al.l (|2007l ) we do have 
a far-UV detection that is offset ~ 0.3" from the primary 
target galaxy and is associated with a faint object in the 
HST optical images. If the object is not a foreground 
contaminant and is at the the same redshift as the tar- 
get galaxy, z = 1.355, the far-UV to optical SED would 
imply no flux decrement shortward of the Lyman limit 
and suggests significant ionizing emission that is signifi- 
cantly displaced from the primary galaxy, similar to some 
of the cand i date Lyman continuum emitters found by 
llwata et al.l (|2009f ). However, we have assumed that the 
souce is simply a faint foreground object. This example 
serves as a caution to LyC studies at higher redshifts, 
where the lines-of-sight are much longer and far more 
subject to contamination (Siana et al. 2007, Vanzella et 
al. 2010, in prep). 

5. DISCUSSION 

We have demonstrated that the average ionizing emis- 
sivity of star-forming galaxies is lower, relative to the 
non- ionizing UV emissivity, at z ~ 1 than at z '-^ 3. Our 
SED fits show that the stellar population parameters of 
our sample do not differ significantly from those of LBGs 
at 2 < z < 3. Therefore, there is no direct evidence 
that the increased LyC emission is due to an intrinsic 
difference in the LyC production from young stars. If 
the stellar populations, and thus LyC productions, are 
the same, that would imply that the increased ionizing 
emission at high redshift is caused by an evolution in the 



escape fraction of ionizing photons due to differences in 
Hi masses and spatial distributions. Of course, it is pos- 
sible that the properties of a subset of the high redshift 
sample may have unique characteristics (top-heavy ini- 
tial mass function, low metallicity) which would increase 
the intrinsic LyC-to-UV ratio, and appear to have a high 
escape fraction. Regardless of whether the escape frac- 
tion is evolving or not, it appears the ionizing emissivity 
is increasing toward higher redshifts. This has important 
implications for models of reionization and the evolution 
of the IGM. 

An evolving Ly C escape fraction has been inferred be- 
fore. Ilnoue et"aLl pOOQ) argue that at z < 2, QSOs can 
provide all of the ionizing background (though the back- 
ground estimates at low redshift are uncertain), so the 
LyC escape fraction from star-forming galaxies at z < 2 
is negligible. However, because of the very low QSO 
space densities at z > 3, they argue that star-forming 
galaxies at very high redshifts must provide the vast ma- 
jority of the ionizing background, which implies higher 
escape fractions at higher redshifts. Of course, this is 
an indirect argument, so it is important to confirm this 
assumption with a direct empirical detection. 

Recent inv estigations of the ionizing backgroun d at 
high redshift ()Bolton et al.ll2005l:lBecker et al.ll2007[ ) note 
that the Hl-photoionization rate is fairly flat between 
2 < z < 4, which implies an incr ease in the ioni zing emis- 
sitivity at higher redshifts (Fauc her-Giguere et al. 2008) . 
However, the total star formation rate density is a ctually 
declining towards higher redshifts iBouwens et al.l ([2C)07i) . 
This seeming contradiction can be resolved if in fact the 
LyC escape fraction is increasing towards higher red- 
shifts. 

Some authors have pointed out the difficulty of ionizing 
the IGM by z ~ 11 (as is sugge sted by the WMAP5 
electron scattering optical depth, iKomatsu et al.l l2009| l 
and maintaining that ionization given the relatively small 
averag e escape fractions, < fesc.r ei >^ 0.1, exhibited at 
z - 3 (|Charvll2008l: iGnedinI [20081) . If the escape fraction 
continues to increase with redshift beyond z = 3, it would 
naturally remedy this discrepancy, elimi nating the need 
to invoke top-heavy IMFs (Charvl [20081 ) or more exotic 
mechanisms. 

It is interesting to compare the evolution in the es- 
cape fraction with predictions from numerical simula- 
tions. This is a difficult problem for simulations, as the 
escape fraction is likely highly sen sitive to the star for - 
mation and feedback presciptions. iGnedin et al.l (|2008l ) 
argue that the escape fraction is highest in high mass 
galaxies. They see no evolution in the escape fraction, 
but they only follow their galaxies from z = 10toz = 4, 
so it is not clear that the y will see an evolving escape 
fraction at lower redshift. [Razoumov fc Sommer-LarsenI 
((20081 get very different results, in that low mass galax- 
ies have the highest escape fractions. Furthermore, they 
see an evolving es cape fraction over all redshifts c overed 
(10 < z < 2.4, Raz oumov fc Sommer-Larsenll20090 . 

Of course, it would be extremely useful to mesure the 
escape fraction at z > 6, near the epoch of reionization. 
However, it is impossible to directly determine if the ion- 
izing emissivity of galaxies continues to increase at z > 3, 
as the increasing IGM opacity makes direct LyC detec- 
tion nearly impossible. In the future, it will be important 
to identify unique characteristics of the LyC emitters at 
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z ^ 3 (at \rest > 1216 A) and determine if galaxies with 
these characteristics are more common during the epoch 
of reionization > 6). 

6. CONCLUSIONS 

We have performed a deep (5 orbits/galaxy) HST 
ACS/SBC far-UV imaging survey of 15 galaxies at z ^ 
1.3 to probe their rest-frame Lyman Continuum (700 A) 
emission. The data achieve a depth of uiABi'^cr) > 28.7 
(AB) in a 1" diameter aperture and are the deepest extra- 
galactic far-UV images of comparable area. Our findings 
are as follows: 

• We do not detect any escaping Lyman Continuum 
from our target galaxies and achieve limits to 
the UV-to-LyC ratio of fi5oo/ hyC > [50,149]. 
Our stacked image gives a 3cr limit of /i5oo//li/C > 
302. After correcting for average IGM opacity (fac- 
tor of 2), and accounting for the intrinsic Lyman 
Break in the SED (factor of 3), these translate to 
individual relative escape fraction limits (3cr) of 

fesc,rel < [0.03,0.21] and fesc,rel < 0.02 in the 

stack. These are the best escape fraction limits 
ever obtained at any redshift. 

• Fits of stellar population models to the rest-frame 
UV to near-IR SEDs of our sample show them to 
have star formation rates, stellar masses, dust ex- 
tinction, and ages similar to LBGs at higher red- 
shift (2 < z < 3), some of which exhibit very high 
escape fractions. We show that the intrinsic Lyman 
continuum flux is extremely difhcult to determine 
as it depends on the star formation history within 
the last < 10 Myr. Therefore, it is impossible to 
infer the absolute escape fraction via SED fits alone 
and other means, such as Ha flux measurements, 
are required. Regardless of the exact value of the 
intrinsic Lyman continuum emissivity, the similar- 
ity in the stellar populations between our sample 
and LBGs suggests that the lack of LyC detections 
at z 3 is not likely due to a relative lack of Ly- 
man continuum production, but rather a difference 
in the Hi mass and its distribution in and around 
galaxies at the two redshifts. 

• In order to properly account for varying opacity of 
the IGM along different lines-of-sight, we perform 
a Monte Carlo simulation to constrain the escape 
fraction with the limits from our data, as well as 
two similar survey s at z 1.3 (Mal kan et al.ll200"3l : 
ISiana et al.ll2007n . At the 95% confidence level, 
there can be no more than 8% of galaxies with very 



high escape fraction {fesc,rei > 0.5). Alternatively, 
if most galaxies have a low, but non-zero, escape 
fraction, it must be less than 0.04. 

• The limits on the LyC-to-UV flu x ratios of th e 
47 galaxies at z ~ 1.3 from .Malkan et al.l (|2003f l: 
ISiana et al.l (|2007f l and this work, are lower than 
any of the LyC detections at z 3. Generally, it 
appears that 8-10% of z 3 galaxies exhibit large 
escape fractions {fesc,rei > 0.5). Our Monte Carlo 
simulations (which account for variations in IGM 
line-of-sight opacity) show that we would have ob- 
served at least one LyC detection 95% of the time 
we perform our experiment. Our lack of any detec- 
tions in such a large sample strongly suggests that 
star-forming galaxies at z ~ 1 have a lower ionizing 
emissivity than galaxies of comparable luminosity 
at z 3. 

• In light of recent claimed detections of escaping 
LyC from areas of low UV surface brightness, we 
have stacked subregions of the galaxies based on 
surface brightness limits (above and below an ob- 
served B-band surface brightness of Sab = 23.04 
mag arcsec"^. We do not detect escaping LyC in 
either sample, with limits of fesc,rei{high) < 0.013 

and fesc,relilow) < 0.03. 

• We do not see evidence of a shadow (on scales 
less than ^ 40 kpc) of the UV background (from 
the Hi in the target galaxies) in a stack of our 
target galaxies. This suggests that either the Hi 
columns are still large at 40 kpc scales or the 
majority of the far-UV background originates in 
the foreground of our sources (eg. z < 1). 

Given the dearth of ionizing emission from z 1 
galaxies, future studies should focus on studying 
properties of z ^ 3 galaxies in detail to determine 
what why some have large escape fractions. In ad- 
dition, it would be useful to identify unique char- 
acteristics of LyC-emitting galaxies at z 3 that 
can be used to identify LyC-emitters at even higher 
redshifts. This would allow for a more accurate es- 
timate of the contribution of star-forming galaxies 
to the high redshift ionizing background and help 
us understand the reionization process in more de- 
tail. 
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